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Abstract Sea-level rise due to both climate change and non-climatic factors threatens coastal settle- 
ments, infrastructure, and ecosystems. Projections of mean global sea-level (GSL) rise provide insufficient 
information to plan adaptive responses; local decisions require local projections that accommodate dif- 
ferent risk tolerances and time frames and that can be linked to storm surge projections. Here we present 
a global set of local sea-level (LSL) projections to inform decisions on timescales ranging from the com- 
ing decades through the 22nd century. We provide complete probability distributions, informed by a 
combination of expert community assessment, expert elicitation, and process modeling. Between the 
years 2000 and 2100, we project a very likely (90% probability) GSL rise of 0.5 -1.2 m under representa- 
tive concentration pathway (RCP) 8.5, 0.4-0.9 m under RCP 4.5, and 0.3-0.8 m under RCP 2.6. Site-to-site 
differences in LSL projections are due to varying non-climatic background uplift or subsidence, oceano- 
graphic effects, and spatially variable responses of the geoid and the lithosphere to shrinking land ice. The 
Antarctic ice sheet (AIS) constitutes a growing share of variance in GSL and LSL projections. In the global 
average and at many locations, it is the dominant source of variance in late 21 st century projections, 
though at some sites oceanographic processes contribute the largest share throughout the century. LSL 
rise dramatically reshapes flood risk, greatly increasing the expected number of "1 -in-10" and "1-in-100" 
year events. 


1. Introduction 

Sea-level rise figures prominently among the consequences of climate change. It impacts settlements 
and ecosystems both through permanent inundation of the lowest-lying areas and by increasing the 
frequency and/or severity of storm surge over a much larger region. In Miami-Dade County, Florida, for 
example, a uniform 90-cm sea-level rise would permanently inundate the residences of about 5% of the 
county's population, about the same fraction currently threatened by the storm tide of a 1 -in-1 00 year 
flood event [Tebaldi et al., 201 2], A 1 -in-1 00 year flood on top of such a sea-level rise would, assuming geo- 
graphically uniform flooding, expose an additional 35% of the population (Climate Central, Surging Seas, 
201 3, retrieved from SurgingSeas.org, updated November 201 3). 

The future rate of mean global sea-level (GSL) rise will be controlled primarily by the thermal expansion 
of ocean water and by mass loss from glaciers, ice caps, and ice sheets [Church et al., 2013]. Changes in 
land water storage, through groundwater depletion and reservoir impoundment, may have influenced 
twentieth-century sea-level change [Gregory et al., 201 3] but are expected to be relatively minor contribu- 
tors compared to other factors in the current century [Church et al., 201 3]. 

Local sea-level (LSL) change can differ significantly from GSL rise [Milne et al., 2009; Stammer et al., 2013], 
so for adaptation planning and risk management, localized assessments are critical. The spatial variability 
of LSL change arises from: (1 ) non-uniform changes in ocean dynamics, heat content, and salinity [Lev- 
ermann et al., 2005; Yin et al., 2009], (2) perturbations in the Earth's gravitational field and crustal height 
(together known as static-equilibrium effects) associated with the redistribution of mass between the 
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cryosphere and the ocean [Kopp et al., 201 0; Mitrovica et at., 201 1 ], (3) glacial isostatic adjustment (GIA) 
[Farrell and Clark, 1 976], and (4) vertical land motion due to tectonics, local groundwater, and hydrocarbon 
withdrawal, and natural sediment compaction and transport [e.g., Miller et al., 201 3]. 

Most past assessments of LSL change have focused on specific regions, such as the Netherlands [Katsman 
etal., 201 1], the U.S. Pacific coast [National Research Council, 201 2], New York City [Horton eta!., 2011; New 
York City Panel on Climate Change, 201 3], and New Jersey [Miller et al., 201 3]. Slangen et al. [201 2], Slangen 
etal. [2014], Perrette et al. [2013], and Church etal. [2013] [AR5] have produced global projections of LSL 
using Coupled Model Intercomparsion Project (CMIP) projections [Taylor et al., 2012] for thermal expan- 
sion and ocean dynamics, along with estimates of net land ice changes, their associated static-equilibrium 
effects and GIA. 

Here we expand upon past efforts to project LSL globally. First, we present a complete probability dis- 
tribution. This is critical for planning purposes; the likely (67% probability) ranges presented in AR5 and 
many other previous efforts provide no information about the highest 1 7% of outcomes, which may well 
be key to risk management. Second, we indicatively extend our projections to 2200, in order to inform 
both decision-making regarding long-term infrastructure investment decisions and their longer term land 
use consequences, and also greenhouse gas mitigation decisions in the context of long-term sea-level rise 
commitments [Levermann et al., 201 3]. Finally, using a Gaussian process model [Kopp, 201 3] of historical 
tide-gauge data [Holgate et al., 201 3], we include probabilistic estimates of local non-climatic factors. 

We first present our framework and projections for selected locations (projections for all tide-gauge loca- 
tions are included in the Supporting Information), then assess the effects of sea-level rise on coastal flood- 
ing risk at these locations. Throughout, we seek to employ transparent assumptions and an easily replica- 
ble methodology that is useful for risk assessment and can be readily updated with new information. 


2. Methods 

LSL projections require the projection and aggregation of the individual components of sea-level change 
[e.g., Milne et al., 2009] at each site of interest. Here, we project three ice sheet components (the Greenland 
Ice Sheet, GIS; the West Antarctic ice sheet, WAIS; and the East Antarctic ice sheet, EAIS); glacier and ice cap 
(GIC) surface mass balance (SMB); global mean thermal expansion and regional ocean steric and ocean 
dynamic effects (which we collectively call oceanographic processes); land water storage; and long-term, 
local, non-climatic sea-level change due to factors such as GIA, sediment compaction, and tectonics. In our 
base case, we allow correlations, derived from the SMB model, between different mountain glaciers but 
otherwise assume that, conditional upon a global radiative forcing pathway, the components are inde- 
pendent of one another. To calculate GSL and LSL probability distributions, we employ 1 0,000 Latin hyper- 
cube samples from time-dependent probability distributions of cumulative sea-level rise contributions 
for each of the individual components. The sources of information used to develop these distributions are 
described below and summarized in Figure 1 . 

We construct separate projections for three representative concentration pathways (RCPs): RCP 2.6, RCP 
4.5, and RCP 8.5 [Meinshausen etal., 2011], which correspond respectively to likely global mean tem- 
perature increases in 2081 -2100 of 1.9-2.3 C, 2.0-3.6 C, and 3.2-5.4 C above 1850-1900 levels [Inter- 
governmental Panel on Climate Change, 201 3]. We do not consider RCP 6.0, as 21 st-century sea-level rise 
projections for this pathway are nearly identical to those for RCP 4.5, and few CMIP Phase 5 (CMIP5) model 
runs for RCP 6.0 extend beyond 2100 [Taylor et al., 201 2], The RCPs do not represent socioeconomic sce- 
narios but can be compared to emissions in no-policy socioeconomic projections such as the Shared 
Socioeconomic Pathways (SSPs) [O'Neill etal., 2014]. Radiative forcing in RCP 6.0 in the second half of 
the century is comparable to that in the lowest emissions SSP (SSP 1 ), while RCP 8.5 is above four of the 
SSPs but below the highest-emission SSP [Riahi, 201 3]. Thus, RCP 8.5 can be viewed as corresponding 
to high-end business-as-usual emissions and RCP 4.5 as a moderate mitigation policy scenario. RCP 2.6 
requires net-negative global emissions in the last quarter of the 21 st century, implying a combination of 
intensive greenhouse gas mitigation and at least modest active carbon dioxide removal. 
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Figure 1. Logical flow of sources of information used in local sea-level projections. GCMs, global climate models; GIC, glaciers and ice 
caps; SMB: surface mass balance. 


2.1. Ice Sheets 

Our projections of 21 st-century changes in mass balance of GIS and the Antarctic ice sheet (AIS) are gener- 
ated by combining the projections of AR5 and the expert elicitation of Bomber and Aspinall [201 3] [BA1 3], 
AR5 is used to characterize median and likely ranges of sea-level change, while BA1 3 is used to calibrate 
the shape of the tails (Supporting Information Figure SI and Table SI). 

AR5 separately assesses AIS and GIS mass balance changes driven by SMB and ice sheet dynamics. For ice 
sheet dynamics, AR5 determined that there was insufficient knowledge to differentiate between RCP 2.6 
and 4.5 (and 8.5 for AIS). Projections of total ice sheet mass loss — given as a likely cumulative sea-level rise 
contribution — are thus partially scenario-independent. BA13 probed more deeply into the tail of ice sheet 
mass loss projections, inquiring into the 5th-95th percentile ranges of GIS, EAIS, and WAIS. However, BA1 3 
does not differentiate between SMB and ice sheet dynamics or between RCPs. 

We reconcile the projections as described in the Supporting Information. For AIS, the reconciled RCP 
8.5 projections (median////ce/y/very likely [90% probability] of 4/-8 to 15/-1 1 to 33 cm) are significantly 
reduced in range relative to BA13 (median///7ce/y/very likely of 13/2 to 41 1-2 to 83 cm); for GIS, the recon- 
ciled projections are almost identical to those based directly on AR5 and have a likely range (8-25 cm) 
close to the very likely range estimated from BA1 3 (9-29 cm) (Supporting Information Table SI). 

Ice sheet mass balance changes do not cause globally uniform sea-level rise. To account for the differing 
patterns of static-equilibrium sea-level rise caused by land ice mass loss, we apply sea-level fingerprints, 
calculated after Mitrovica et al. [201 1] (Supporting Information Figure S2). These fingerprints assume mass 
loss from each ice sheet is uniform; in most regions, the error introduced by this assumption is minimal 
[Mitrovica et al., 201 1 ]. 

2.2. Glacier and Ice Caps 

For each RCP, we generate mass balance projections for 1 7 different source regions of glaciers and ice caps 
(described in the Supporting Information). For each source region, we employ a multivariate f-distribution 
of ice mass change with a mean and covariance estimated from the process model results of Marzeion 
et al. [201 2]. Each source region has a distinct static-equilibrium sea-level fingerprint, calculated in the 
same fashion as for ice sheet mass loss (Supporting Information Figure S2). 

The projections based on Marzeion etal. [2012] are modestly narrower and have a slightly higher median 
than those of AR5: a likely range of 9- 1 5 cm from non-Antarctic glaciers by 21 00 for RCP 2.6 (vs. 4- 1 6 cm 
for AR5) and 14-21 cm for RCP 8.5 (vs. 9-23 cm for AR5). We opt for the Marzeion etal. [2012] projections 
because of the availability of disaggregated output representing projections based on a suite of global 
climate models (GCMs) for each source region. 
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2.3. Oceanographic Processes 

Projections of changes in GSL due to thermal expansion and in LSL due to regional steric and dynamic 
effects are based upon the CMIP5 GCMs. In particular, we employ a t-distribution with the mean and 
covariance of a multi-model ensemble constructed from the CMIP5 archive (Supporting Information 
Figures S3 and S4, Table S2). Values used are 1 9-year running averages. For each model, we use a single 
realization. The sea-level change at each tide-gauge location is assumed to be represented by the nearest 
ocean grid cell value of each GCM. 

The horizontal resolution of the CMIP5 ocean models is ~ 1 degree. In these coarse-resolution models, 
sea level at the coast may differ from the open ocean due to local biases driven by unresolved processes 
(e.g., coastal currents) and bathymetry [Holt et al., 2009] or via the influence of small-scale processes (e.g., 
eddies) on larger-scale steric and dynamic changes [Penduff et al., 2010, 2011]. Although there is some 
evidence that climate-forced trends in sea level are not sensitive to resolution [Penduff et al., 201 1; Suzuki 
et al., 2005], higher-resolution coastal modeling is required to determine whether the probabilities esti- 
mated at the GCM grid scale are significantly changed by sub-grid processes. 

GCM projections exhibit a range of late nineteenth-century sea-level behavior largely attributable to 
model drift. Uncorrected GCM-based estimates of the rate of mean global sea-level change from 1861 to 
1 900 range from -0.4 to +1 .1 mm/yr. To correct for global-mean model drift, we apply a linear correction 
term to each model. The linear correction adjusts the rate of GSL rise over 1861-1900 to match a rate of 
thermal expansion estimated by removing the multi-model average of GIC mass loss from Marzeion et al. 
[2012] from the GSL curve of Church and White [201 1], After correction, the rate of thermal expansion over 
1861-1 900 is 0.3 ± 0.9 (2rr) mm/yr (Supporting Information Figure S3). 

Consistent with AR5's judgment that the 5th-95th percentile of CMIP5 output represents a likely (67% 
probability) range for global mean thermal expansion, we multiply the standard deviation of the 
f-distribution for oceanographic processes by 1 .7. 

2.4. Land Water Storage 

Following the approach of Rahmstorf et al. [201 2], we estimate GSL change due to changes in water stor- 
age on land based upon the relationship between such changes and population (Supporting Information 
Figure S5). For changes in reservoir storage, we use the historical cumulative impoundment estimate of 
Chao et al. [2008]. We assume that reservoir construction is a sigmoidal function of population: 


where / is impoundment expressed in mm equivalent sea level (esl), P(t) is world population as a function 
of time, and the remaining variables are constants. The results imply a maximum additional impound- 
ment of 6 mm (esl) on top of the current 30 mm; based on the discrepancy between the "nominal" and 
"actual" impoundment estimated by Chao et al. [2008], we conservatively allow a 2a error in this estimate 
of ±50%. 

For the rate of groundwater depletion, we fit the estimates of Wada et al. [2012] and Konikow [2011] as 
linear functions of population, forced through the origin. The estimate of Wada et al. [201 2] is based 
on fluxes estimated from a global hydrological model of groundwater recharge and a global database 
of groundwater abstraction, while that of Konikow [201 1] uses a range of approaches depending on 
the data available for each aquifer. We take the mean and standard deviation of the two slopes esti- 
mated (0.06 ± 0.02 mm/yr/billion people) and allow an additional 2a error of ± 50 %, a level based 
upon the errors estimated by the authors of the two impoundment studies. In our main calcula- 
tion, we do not include the water resource assessment model-based estimate of Pokhrel et al. [201 2], 
which is about a factor of three higher than the other two estimates; we include this estimate in a 
sensitivity case. 

We employ population projections derived from United Nations Department of Economics and Social 
Affairs [2014]. We treat population as distributed following a triangular distribution, with the median, 
minimum, and maximum values corresponding to the middle, low, and high U.N. scenarios (10.9, 6.8, and 
16.6 billion people in 2100, respectively). For scenarios in which population declines, we allow 


I = ax erf((P(f) - b) /c) + l 0 


( 1 ) 
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some reduction in impoundment, but do not allow impoundment to decrease below its year 2000 
level. 

2.5. Glacial Isostatic Adjustment, Tectonics, and Other Non-Climatic Local Effects 

GIA, tectonics, and other non-climatic local effects that can be approximated as linear trends over the 
twentieth century are assumed to continue unchanged in the 21st and 22nd centuries. This is a good 
assumption for GIA, but imperfect for other processes. Tectonic processes can operate unsteadily, and a 
linear trend estimated from the historical record may be inaccurate. LSL rise related to fluid withdrawal is 
subject to engineering, resource depletion, market factors and policy controls, and might either increase 
or decrease in the future relative to historical levels. In addition, the trend estimates can encompass slow 
ocean dynamic changes that are close to constant over the historical record but could change in the 
future. Nonetheless, fora global analysis, assuming the continuation of observed historical changes offers 
the best currently feasible approach. 

We estimate historical rates using a spatiotemporal Gaussian process model akin to that employed by 
Kopp [2013]. Sea level as recorded in the tide-gauge records (Permanent Service for Mean Sea Level, 
Tide-gauge data, retrieved from http://www.psmsl.org/data/obtaining/, accessed January 2014) is rep- 
resented as the sum of three Gaussian processes: (1 ) a globally uniform process, (2) a regionally varying, 
temporally linear process, and (3) a regionally varying, temporally autocorrelated non-linear process. We 
allow for spatial non-stationarity in the Gaussian process prior by optimizing the hyperparameters sepa- 
rately for each of 1 5 regions (Supporting Information Table S4 and Figure S6). The posterior estimate of 
the second (linear) process at each site is used for forward projections. Mathematical details are provided 
in the Supporting Information. 

2.6. Post-2100 Projections 

Indicative post-2100 projections are developed according to the methods described in the previous 
sections. For ice sheet mass balance, we continue the constant 21st century acceleration. For non-climatic 
factors that are approximated as linear in the 21st century, we continue the constant 21st century rate. For 
land water storage, we extend the population projections using the 22nd century growth rates of United 
Nations Department of Economics and Social Affairs [2004] and use the same relationships of impound- 
ment and groundwater depletion to population as in the 21st century (Supporting Information Figure 
S5). The number of GCM-based model results for GIC and oceanographic processes drops significantly 
beyond 2100 (Supporting Information Table S2), leading in these terms to a modest discontinuity and a 
reduction in variance in these terms at the start of the 22nd century (Supporting Information Figure S7). 
Acknowledging the limitations of these assumptions, we present post-21 00 projections in tables rounded 
to the nearest decimeter. 

2.7. Flood Probabilities 

To examine the implications of our projections for coastal flooding, we combine Latin hypercube samples 
from the sea-level distribution for an illustrative subset of sites with maximum-likelihood generalized 
Pareto distributions (GPDs) estimated from observed storm tides after Tebaldi et al. [2012], updated to 
use the full historic record of hourly water levels available at each location. Flourly data for non-U. S. sites 
are from the University of Flawaii Sea Level Center (retrieved from uhslc.soest.hawaii.edu, May 2014). 

The estimated GPDs do not take into account any future changes in storm frequency, intensity, or track 
[e.g., Knutson eta!., 2010], so projected future flood probabilities should be viewed primarily as an 
illustration. 

Using the maximum-likelihood GPDs, we compute return levels corresponding to a set of representative 
return periods (e.g., the 1 -in-10 or 1 -in-1 00 year flood events). For each decade of each realization of LSL 
change, we then add the projected sea-level change and re-estimate a GPD. The result for each realization 
is a trajectory of probabilities over time for each of the original return levels. For example, for the 10-year 
event, the initial probability at 2000 is 1 0% per year and increases over time as sea-level rises. Cumula- 
tively summing each decade's expectation through the century, we compute the expected number of the 
original events by 2100. Under stationary sea levels, there would be one expected 1 -in-1 00 year event and 
ten expected 1 -in-10 year events between 2001 and 2100. 
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Figure 2. Projections of cumulative contributions of (a) the Greenland ice sheet, (b) the Antarctic ice sheet, (c) thermal expansion, 
and (d) glaciers to sea-level rise in RCP 8.5. Heavy = median, light = 67% range, dashed = 5th-95th percentile; dotted = 0.5th -99.5th 
percentiles. 


3. Sea-Level Projections 

3.1. Mean Global Sea-Level Projections 

The cumulative 21st century GSL contribution of each component is shown in Figure 2 (for RCP 8.5) and in 
Table 1 and Supporting Information Figure S7 (for all RCPs). In the 21 st century, thermal expansion and GIC 
provide the largest contributions to the median outcome and have narrower uncertainty ranges than the 
ice sheet contributions. AIS has the broadest uncertainty range, extending from a small negative contribu- 
tion to sea level (presumably due to warming-induced increased snow accumulation) to a large positive 
contribution (requiring a substantial and/or widespread dynamic change). 

Adding samples from the component distributions together indicates a likely GSL rise (Figure 3 and 
Table 1) in RCP 8.5 of 0.6- 1.0 m by 2100, with a very likely range of 0.5-1 .2 m and a virtually certain (99% 
probability) range of 0.4- 1 .8 m. The right-skewed "fat tail" of the projections arises from the ice sheet 
components. Even in the low-emissions RCP 2.6 pathway, sea-level rise by 2100 very likely exceeds the 
32 cm that would be projected from a simple linear continuation of the 1993-2009 rate [Church and 
White, 2011]. 

Through the middle of the current century, GSL rise is nearly indistinguishable between the three forcing 
pathways (Figure 3 and Table 1). Only in the second half of the century do differences of >6 cm begin to 
emerge in either the median or the tails of the projections. By 2100, median projections reach 0.8 m for 
RCP 8.5, 0.6 m for RCP 4.5 and 0.5 m for RCP 2.6. By 2200, upper tail outcomes are clearly higher in the 
high-forcing pathway, yet there remains significant overlap in the ranges of all three pathways, with likely 
GSL rise by 2200 of 1. 3-2.8 m in RCP 8.5 and 0.5- 1.6 m in RCP 2.6. The overlap between RCPs is due in 
significant part to the large and scenario-independent uncertainty of AIS dynamics, even as the thermal 
expansion, GIC and, to a lesser extent, GIS contributions begin to differentiate (Supporting Information 
Figure S7). 
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Table 1 . GSL Projections 3 



RCP 8.5 

RCP 4.5 

RCP 2.6 

cm 

50 

17-83 

5-95 

0.5-99.5 

99.9 

50 

17-83 

5-95 

0.5-99.5 

99.9 

50 

17-83 

5-95 

0.5-99.5 

99.9 

2100— 

Components 














GIC 

18 

14-21 

11-24 

7-29 

<30 

13 

10-17 

7-19 

3-23 

<25 

12 

9-15 

7-17 

3-20 

<25 

GIS 

14 

8-25 

5-39 

3-70 

<95 

9 

4-15 

2-23 

0-40 

<55 

6 

4-12 

3-17 

2-31 

<45 

AIS 

4 

-8 to 15 

-11 to 33 

-14 to 91 

<155 

5 

-5 to 16 

-9 to 33 

- 1 1 to 88 

<150 

6 

-4 to 17 

-8 to 35 

-10 to 93 

<155 

TE 

37 

28-46 

22-52 

12-62 

<65 

26 

18-34 

13-40 

4-48 

<55 

19 

13-26 

8-31 

1-38 

<40 

LWS 

5 

3-7 

2-8 

- 0 to 1 1 

<11 

5 

3-7 

2-8 

-Oto 11 

<11 

5 

3-7 

2-8 

-Oto 11 

<11 

Total 

79 

62-100 

52-121 

.39- 1 76 

<245 

59 

45-77 

36-93 

24-147 

<215 

50 

37-65 

29-82 

19-141 

<210 


Projections by year 


2030 

14 

12-17 

11-18 

8-21 

<25 

14 

12-16 

10-18 

8-20 

<20 

14 

12-16 

10-18 

8-20 

<20 

2050 

29 

24-34 

21-38 

16-49 

<60 

26 

21-31 

18-35 

14-44 

<55 

25 

21-29 

18-33 

14-43 

<55 

2100 

79 

62-100 

52-121 

39-176 

<245 

59 

45-77 

36-93 

24-147 

<215 

50 

37-65 

29-82 

19-141 

<210 

2150 

130 

100-180 

80-230 

60-370 

<540 

90 

60-130 

40- 1 70 

20-310 

<480 

70 

50-110 

30-150 

20-290 

<460 

2200 

200 

130-280 

100-370 

60-630 

<950 

130 

70-200 

40-270 

10-520 

<830 

100 

50-160 

30-240 

10-500 

<810 


Other projections for 2100 


AR5 

73 

53-97 


52 

35-70 

43 

28-60 


H14 


70-120 

50- 1 50 




40-60 

25-70 

J12 

110 


81-165 

75 

52-110 

57 


36-83 

SI 2 




90 

64-121 

75 


52-96 


a TE: Thermal expansion, LWS: Land water storage, HI 4: Horton etal. [2014], J 1 2: Jevrejeva et al. [2012], SI 2: Schaeffer et ai. [2012]. 
All values are cm above 2000 CE baseline except for AR5, which is above a 1986-2005 baseline. 


The importance of different compo- 
nents to the GSL uncertainty varies with 
time. While in 2020 about two-thirds 
of the total variance in GSL is due to 
uncertainty in projections of thermal 
expansion, by 2050 in RCP 8.5 changes 
in ice sheet volume are responsible 
for more than half the variance and 
changes in thermal expansion for only 
about one-third. By 2100, AIS alone 
is responsible for half the variance, 
with an additional 30% due to GIS 
uncertainty and only 1 5% due to uncer- 
tainty in thermal expansion (Figure 4). 
Because the uncertainty in AIS mass 
loss is largely scenario-independent, 
its dominant contribution to variance 
holds across RCPs; indeed, it is even 
more dominant in lower emissions pathways where the contributions from other sources are smaller and 
more strongly constrained (Supporting Information Figure S8). 

3.2. Comparison With Other Global Projections 

By construction, our likely projections of GSL in 2100 are close to those of AR5 (Table 1), though differ 
slightly (e.g., in RCP 8.5 in 2100, 0.6- 1.0 m vs. AR5's 0.5 -1.0 m) due to: (1) the drift correction to a possibly 
non-zero (0.3 ± 0.9 mm/yr) background thermal expansion, (2) the use of Marzeion et al. [2012] for GIC, 
and (3) the use of a year 2000 as opposed to 1985-2005 baseline. AR5 projections of GSL rise are lower 



2000 2050 2100 2150 2200 


Figure 3. Projections of GSL rise for the three RCPs. Heavy = median, 
dashed = 5th-95th percentile, dotted = 0.5th -99.5th percentiles. 
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Figure 4. Sources of variance in raw (a, c) and fractional terms (b, d), globally (a — b) and at New York City (c-d) in RCP 8.5. AIS: 
Antarctic ice sheet, GIS: Greenland ice sheet, TE: thermal expansion. Ocean: oceanographic processes, GIC: glaciers and ice caps, LWS: 
land water storage, Bkgd: local background effects. 


than those from other sources, such as semi-empirical models [Rahmstorf, 2007; Schaeffer et al., 2012; 
Vermeer and Rahmstorf, 2009] and expert surveys [Horton eta!., 2014]. However, AR5 only projects likely 
ranges; higher magnitudes of ice loss are implied if less likely outcomes are considered [Little et al., 201 3a], 

By using plausible information to complement the AR5 analysis, we project a very likely GSL rise in 21 00 
of 0.4-0.9 m for RCP 4.5, which compares to the 90% probability semi-empirical projections of 0.5 -1.1 m 
[Jevrejeva etai, 2012] and 0.6 -1.2 m [Schaeffer et al., 201 2], The widths of the semi-empirical very likely 
ranges are similar to those of our projections, with the entire distribution shifted to higher values. The 
95th percentiles of these two semi-empirical projections resemble the 98th and 99th percentiles of our 
projection, respectively. 

Horton etal. [2014] conducted a survey of 90 experts with a substantial published record in sea-level 
research. Their survey found likely /very likely sea-level rise by 2100 of 0.7-1. 2/0.5 -1.5 m under RCP 8.5 and 
0.4-0.6/0.3-0.7 m under RCP 2.6. Our projections for RCP 2.6 are similar to those of the surveyed experts, 
with a slightly fatter upper tail, while the experts' responses for RCP 8.5 are considerably fatter-tailed than 
our projections. The surveyed experts' 83rd and 95th percentiles correspond to our 95th and 99th per- 
centiles, respectively. Although Horton etal. [2014] did not probe the reasons why their surveyed experts 
differed from AR5, we suggest it may be related to expectations about the behavior of Antarctica. As 
noted previously, high-end estimates of Antarctic mass loss from the expert elicitation of BA13 are higher 
than would be expected from the likely range of AR5 projections; our reconciled ice sheet projections 
significantly lower this contribution. (See also the sensitivity tests in section 4 for comparison.) 

Our 99.9th percentile estimate for 2100 under RCP 8.5, 2.5 m, is consistent with other estimates of the 
maximum physically possible rate of sea-level rise in the 21 st century [e.g., Miller et al., 201 3]. It is higher 
than the 2.0 m estimate of Pfeffer et al. [2008], which corresponds to our 99.7th percentile. Comparing 
the individual contribution to Pfeffer et a l.'s high-end projection shows that their projected GIC mass loss 
(55 cm) exceeds our highest projected value (32 cm), while their projections of GIS and AIS mass loss 
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(a) Median projection: RCP 8.5 

GSL = 0.79 m 
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• • 0.45 

[ J 0.40 


(53 and 62 cm, respectively) correspond to our 98th and 99th percentiles. Their projection of thermal 
expansion (30 cm) includes no uncertainty and corresponds to our 22nd percentile [cf., Sriver et a!., 201 2]. 
They do not include changes in sea-level rise resulting from changes in land water storage. As noted pre- 
viously, the tail of the sea-level rise projections is dominated by the uncertainty in AIS mass loss, which is 
lower in Pfefferetal. [2008] than in either our projections or BA13. 

3.3. Local Sea-Level Projections: Patterns 

Figure 5 displays the median LSL projections for RCP 8.5 in 2100 and the projection uncertainty, as 
reflected by the difference between the 17th and 83rd percentile levels. In Figure 6, we illustrate the 
relationship between LSL and GSL using three indices: (1 ) the median value of R, which we define as the 
ratio of LSL change driven by land ice and oceanographic components to GSL change driven by those 
same components (Figure 6a), (2) the uncertainty in R, reflected in Figure 6b by the difference between 
its 1 7th and 83rd percentile levels, and (3) the magnitude and uncertainty of background, non-climatic 
LSL change (Figure 6c). For sites where R is close to 1 and exhibits little uncertainty, GSL projections with 
adjustment for local land motion provide a reasonable estimate of LSL; for other sites, more detailed 
projections, such as those in this article, are necessary. 

The median value of R (Figure 6a) is within 5 % of unity at about a quarter of tide-gauge sites, with higher 
values in much of Oceania, the Indian Ocean, and southern Africa resulting from the static-equilibrium 
effects of land ice mass loss. R generally declines toward higher latitudes due to static-equilibrium effects, 
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Background rate (mm/y) 




Figure 6. (a) Median ratio R of climatically driven LSL change to climatically driven GSL exchange (i.e., excluding land water storage 
and local land motion) in RCP 8.5 in 2100; (b) width of the 17th -83rd percentile range of R ; (c) mean estimates of background rate of 
sea-level rise due to GIA, tectonics, and other local factors (mm/yr). Open circles in bottom indicate sites where 2 n range spans zero. 
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but with values elevated in northeastern North America and to a lesser extent the North and Baltic Seas 
by oceanographic processes. This pattern — with sea-level rise dampened near land ice and enhanced far 
from it and in the northwestern North Atlantic — resembles that found by previous studies [Kopp et al., 
2010; Perrette et al., 2013; Slangen et al., 2012, 2014], Uncertainty in R (Figure 6b) is also relatively small 
(likely range width of < 30 %) in the inhabited southern hemisphere and low-latitude northern hemi- 
sphere, with the range increasing northwards due to both the sensitivity of static-equilibrium effects to 
the particular distribution of shrinking land ice reservoirs and — especially in northeastern North America, 
the Baltic Sea, and the Russian Arctic — uncertainty in oceanographic processes (Supporting Information 
Figure S4). 

Added on top of the climatically driven factors reflected in Rare the global effects of land water storage 
(not shown in Figure 6) and the effects of local land motion (Figure 6c). Moderately high rates of land sub- 
sidence can be associated with GIA, as in the northeastern United States (e.g., 1 .3 ± 0.2 mm/yr at New 
York City), while more extreme rates generally include contributions from fluid withdrawal, delta pro- 
cesses, and/or tectonics. Subsidence driven by fluid withdrawal and delta processes is high at sites such as 
Bangkok, Thailand (background rate of 1 1 .9 ± 1 .1 mm/yr at the Fort Phracula Chomklao tide gauge), Grand 
Isle, Louisiana (7.2 ± 0.5 mm/yr), Manila, the Philippines (background rate of 4.9 ± 0.6 mm/yr), and Kolkata, 
India (5.1 ± 1 .0 mm/yr). Episodic tectonic factors play an important role in both subsidence and uplift 
in Japan, where average background rates can range from - 5.2 ± 0.7 mm/yr at Onahama to 1 8.0 ± 1 .6 
mm/yr at Toba. At high latitudes, GIA-related uplift gives rise to high background rates of sea-level fall, as 
can be seen in places like Juneau, Alaska, (- 14.9 ±0.5 mm/yr), and Ratan, Sweden (- 9.3 ±0.2 mm/yr). 
While some previous global projections have used physical models to incorporate GIA [e.g., Slangen et al., 
201 2, 2014], the current projections are to our knowledge the first to employ observationally based rates. 

3.4. Local Sea-Level Projections: Examples 

To illustrate the importance of local factors for sea-level rise projections, we consider several sites along 
the coasts of the United States where different factors dominate LSL change (Tables 2 and 3). Vv/hile we 
focus on projections for RCP 8.5 as a way of highlighting the differences between GSL and LSL, similar 
considerations apply to other RCPs, which are shown in the tables. 

New York City experiences greater-than-global sea-level rise under almost all plausible projections, with a 
likely range of 0.7-1 .3 m by 2100 under RCP 8.5. Three factors enhance sea-level rise at New York. First, due 
to its location on the subsiding peripheral bulge of the former Laurentide Ice Sheet, the site experiences 
GIA-related sea-level rise of 1.3 ±0.2 mm/yr. Second, the rotational effects of WAIS mass loss increase 
the region's sea-level response to WAIS mass loss by about 20% [Mitrovica et al., 2009]. Third, as noted 
in earlier papers [Kopp et al., 201 0; Yin et al., 2009; Yin and Goddard, 201 3], changes in the Gulf Stream may 
result in dynamic sea-level rise in the mid-Atlantic United States. This enhancement can be seen by exam- 
ining the difference between oceanographic sea-level rise at New York and the global average, which 
has a median of 14cm and a likely range of -6 to +35 cm. These three effects are partially counteracted 
by the — 55 % reduction in the sea-level response due to GIS mass loss, associated with the gravitation- 
ally induced migration of water away from of this relatively proximal ice mass. Indeed, the climatic fac- 
tors that amplify and reduce LSL rise relative to GSL rise are nearly balanced in the median projection 
(/?= 1.03, with a likely range of 0.73-1.30), with GIA effects pushing local rise to levels that exceed the 
global rise. 

Sewell's Point in Norfolk, VA, is projected to experience higher-than-global mean sea-level rise due to 
the same factors as New York City: subsidence due to GIA, enhanced influence of WAIS mass loss, and 
exposure to changes in the Gulf Stream. Being located farther south along the U.S. East Coast, Norfolk 
experiences somewhat smaller ocean dynamic changes (median and likely ocean dynamic sea-level rise 
increment of 9 cm and -8 to 26 cm) but greater sea-level rise due to GIS mass loss (experiencing about 
~ 45 % less sea-level rise than the global mean). Its R value (1.00, likely range of 0.75-1.22) is similar to 
New York City. However, whereas New York City sits upon bedrock, Norfolk is located on the soft sedi- 
ments of the Coastal Plain [Miller et al., 201 3]. As a consequence, it is exposed to sea-level rise due to both 
natural sediment compaction and compaction caused by groundwater withdrawal, which increases the 
background non-climatic rate of sea-level rise to 2.5 ± 0.3 mm/yr. Accordingly, the likely range of LSL rise 
for RCP 8.5 in 21 00 is 0.8- 1.3 m. 
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Sea-level rise at Key West, Florida, is closer to the global mean, with a likely range in RCP 8.5 by 21 00 of 
0.6-1 .1 m (median R = 1.00, likely range of 0.83-1.15, background rise of 0.5 ± 0.4 mm/yr). By contrast, the 
deltaic western Gulf of Mexico coastline experiences some of the fastest rates of sea-level rise in the world 
as a result of groundwater withdrawal and hydrocarbon production [Kolker et a!., 2011; White and Trem- 
blay, 1 995]. At Galveston, Texas, a background subsidence rate of 4.6 ± 0.3 mm/yr drives a likely range of 
sea-level rise by 21 00 in RCP 8.5 of 1 .0- 1 .5 m. Because the uncertainty in subsidence rate is small relative 
to other sources of uncertainty, this causes a shift in the range rather than a broadening of overall uncer- 
tainty, as occurs at New York City (reflected in a likely R of 0.78-1.13, which is narrower than at New York 
City). 

The Pacific Coast of the contiguous United States is subject to considerable short length-scale sea-level 
rise variability due to tectonics, as can be seen by comparing the background non-climatic rate of 
sea-level rise at Los Angeles (- 1 .1 ± 0.3 mm/yr) and nearby Santa Monica (- 0.6 ± 0.3 mm/yr). In general, 
sea-level rise on this coast is close to the global average, with a likely range in 2100 under RCP 8.5 at 
San Francisco of 0.6-1 .0 m (median R = 0.96, likely 0.84-1 .08, background rate of - 0.1 ± 0.2 mm/yr). 

The slightly lower-than-global projection is a result of smaller Greenland and GIC contributions due to 
proximity to these land ice reservoirs, though counterbalanced by enhanced sea-level rise from AIS mass 
loss. Ocean dynamic factors are projected to play a minimal role. 

Farther north, the proximity of historic and modern glaciers controls LSL projections. At Juneau, pre- 
dicted sea-level rise is dominated by a glacio-isostatic sea-level fall of 14.9 ±0.5 mm/yr, interpreted as 
resulting primarily from the ongoing response to post-Little Ice Age glacial mass loss, with a secondary 
contribution from post-Last Glacial Maximum GIA [Larsen eta!., 2005]. Moreover, shrinking glaciers in 
Alaska and western Canada cause about 2.4 mm of LSL fall at Juneau for every mm of global sea-level 
rise, which reduces the overall magnitude of sea-level rise caused by projected glacial mass loss (median 
R = 0.71 , likely 0.59-0.83). As a consequence, under RCP 8.5 Juneau is likely to experience a sea level fall of 
0.7-1. 1 m by 2100. 

Hawaii and other central Pacific islands experience significantly greater-than-average sea-level rise result- 
ing from land ice mass loss (20% enhancement for GIS, EAIS, and the median combination of shrinking 
glaciers, and 30% for WAIS, giving rise to median R= 1.13 and likely 0.98 -1.26). The likely range of sea-level 
rise at Honolulu, Hawaii, is slightly higher than the global mean (0.6- 1.1 m in 2100 under RCP 8.5, with a 
background rate of - 0.2 ± 0.4 mm/yr). The amplification relative to the global mean is more apparent in 
the tail of the projections, where ice sheet mass loss contributions constitute a larger proportion of the 
sea-level rise. As a consequence, the tail of sea-level rise is fatter at Hawaii than globally, with a 95th per- 
centile in RCP 8.5 of 1 .4 m (compared to GSL of 1 .2 m) and a 99.5th percentile of 2.1 m (compared to GSL 
of 1.8 m). 

A similar range of behaviors is seen outside the United States. At Cuxhaven, on the German North Sea 
coast, a slightly higher-than-global likely range of 0.6-1 .1 m arises from a background subsidence rate of 
1 .0 ± 0.2 mm/yr. Because of its relative proximity to Greenland, Cuxhaven is less exposed to climatically 
driven sea-level rise than average (median R = 0.89, likely 0.62-1.15); unlike sites in eastern North Amer- 
ica that are similarly close to Greenland, it does not experience a countervailing oceanographic sea-level 
rise. The city of Stockholm, Sweden, like Juneau, is experiencing a strong GIA-related uplift of - 5.0 ±0.1 
mm/yr, leading to a likely sea-level rise of -0.4 to +0.8 m. Being farther from a large, actively shrinking 
glacier, however, Stockholm is in the median more exposed than Juneau to climatically driven sea-level 
change (median R = 0.83, likely 0.41 - 1 .20). 

Like Honolulu, the town of Kushimoto, in Wakayama Prefacture, Japan, is in the far-field of the major ice 
sheets and most major glaciers. It is also exposed to a likely ocean dynamic sea-level rise of -5 to +1 8 cm 
in 2100. Together, these factors lead to a median R= 1.14, likely 0.98-1 .28. Kushimoto also is experiencing 
tectonic subsidence, leading to a likely sea-level rise in 2100 of 0.8- 1.3 m. 

The city of Valparaiso, on the Chilean Pacific coast, is experiencing tectonic uplift of 2.5 ± 0.8 and exposed 
to a likely ocean dynamic sea-level fall of -2 to 9 cm. Although it experiences about 30% less-than-global 
sea-level rise due to WAIS mass loss, it experiences a larger-than-average response to GIS, EAIS, and most 
glaciers; accordingly its overall sensitivity to sea-level rise is close to the global average (median R = 0.99, 
likely 0.90-1 .08). All these factors together yield a likely 2100 sea-level rise 0.4-0.8 m. 
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4. Variance and Sensitivity Assessment 

As shown in the previous section, LSL rise is controlled by different factors — both climatic and 
non-climatic — at different locations and intervals over the next two centuries. The analysis also reveals 
that the adopted risk tolerance (choice of exceedance probability) influences the importance of different 
components. Median outcomes will vary regionally, driven strongly by varying levels of subsidence and, in 
certain regions, oceanographic processes. High-end (low-probability) outcomes are driven, globally and in 
most locations, by uncertainty in the ice sheet contribution, with the Antarctic signal becoming dominant 
in the highest end of the tail, particularly later in the century (Figure 4 and Supporting Information Figures 
S8 and S9). This contribution varies less by location. 

To test the robustness of our results, we examine three alternate assumptions regarding ice sheet mass 
loss and two alternative assumptions regarding the robustness of GCM projections (Supporting 
Information Tables SI and S3): 

1 . AR: using a lognormal fit to the AR5 median and likely ranges of ice sheet mass balance (GIS almost 
unchanged from reconciled projections; for AIS, very likely range of -15 to 23 cm in RCP 8.5 by 2100) 

2. BA: using a lognormal fit to the BA1 3 median and very likely projections of ice sheet mass balance 
(GIS: median 14 cm and very likely 9-29 cm; AIS: median 14 cm, very likely -2 to 83 cm) 

3. Alt. Corr.: assuming positive correlations of 0.7 between WAIS and GIS and a negative correlation of 
-0.2 between EAIS and the other two ice sheets, following the main projections of Bamber and 
Aspinall [2013] 

4. High GCM Confidence: assuming the very likely ranges estimated by the GCMs for oceanographic 
changes are very likely rather than likely ranges 

5. Reduced degrees of freedom (DOF): Assuming the GCMs collectively provide only six independent 
estimates of GIC and oceanographic change, due to non-independence of models 

6. Higher groundwater depletion (GWD): The ratio of groundwater depletion to population is treated as 
a triangular distribution, with the minimum, median, and maximum estimated respectively from 
/Con//cow [201 1], Wada etal. [2012] and Pokhrel et al. [2012]. 

At a global level and at most locations, the two alternative characterizations of ice sheet mass changes 
have the largest effects, with the median sea-level rise under RCP 8.5 in 2100 varying between 79cm 
under default assumptions and case AR and 91 cm under case BA. The effect is larger in the tails, with 
99.5th percentile projections of 140 cm under AR, 176 cm under default assumptions, 187 cm under Alt. 
Corr., and 300 cm under BA. Varying the confidence in GCMs, by contrast, has little global effect. Although 
the per-capita rate of groundwater depletion estimated from Pokhrel et al. [201 2] is about three times that 
of the Wada etal. [2012], the overall effect of the Higher GWD assumption is small due to the magnitude 
of other uncertainties; this case experiences 3 cm extra GSL rise at the 5th percentile, 4 cm at the median, 
and 6 cm at the 99.5th percentile. 

While all LSL projections are sensitive to assumptions about ice sheet behavior, some are sensitive to 
assumptions about confidence in GCM output. Due to the wide range of projections in the CMIP5 ensem- 
ble at New York City [Yin, 2012], the 99.5th percentile projections are 21 2 cm under default assumptions, 
205 cm under High GCM Confidence, and 232 cm under Reduced DOF. Even at New York City, however, 
GCM uncertainty remains secondary to ice sheet uncertainty; the 99.5th percentile is 1 78 cm under AR, 
212 cm under Alt. Corr., and 359 cm under BA. Moreover, the significance of GCM uncertainty can be quite 
small: at the sites discussed above, the difference between the 99.5th percentiles of the High GCM Confi- 
dence and Reduced DOF cases under RCP 8.5 in 2100 is 27 cm at New York, 19 cm at Sewell's Point, 12 cm 
at Cuxhaven, 9 cm at Galveston, and 6 cm or less at Honolulu, Juneau, Key West, Kushimoto, San Fran- 
cisco, and Valparaiso. A large difference (55 cm) at Stockholm may reflect differences between GCMs in 
the representation of the semi-closed Baltic Sea. 

These sensitivity analyses are not exhaustive. There remains a need for improved ice sheet models to allow 
robust projections of the ice sheet component without heavy reliance upon expert elicitation. However, 
the development of such models is hindered by the limited consensus on the magnitude of positive 
and negative feedbacks on ice loss, such as those involving (a) temperature and snow albedo [Picard 
et al., 201 2], (b) forest fires and snow albedo [Keegan et al., 2014], (c) snowfall and ice sheet discharge 
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[Winkelmann et al., 201 2], (d) grounding line retreat [Joughin et at., 201 4; Rignot et at., 201 4; Schoof, 2007], 
(e) static-equilibrium sea-level and grounding line retreat [Gomez et at., 2010, 2012, 2013], (f) meltwa- 
ter, ocean temperature, sea ice, and snowfall [Bintanja et at., 201 3], and (g) ice-cliff collapse [Bassis and 
Walker, 201 2; Pollard and DeConto, 2013]. The wide range of projections and underlying uncertainties in 
continental-scale model projections pose challenges for interpreting the likelihood of their results [Bind- 
schadler et al., 2013]. It is possible, however, that incomplete information could be better integrated in a 
probabilistic framework [Little et al., 201 3a, 201 3b]. 

Furthermore, structural errors in models of other sea level components remain probable. These errors 
(e.g., a systematic bias caused by a missing process and/or feedback) may have a large impact on tails. 
Here, we do not attempt to perform a systematic analysis. However, we believe that this framework may 
be used to effectively allow for these possibilities to be considered. The subjective judgment applied in 
formulating these distributions is explicit and may be revisited overtime. 

5. Implications for Coastal Flooding 

Since our projections provide full probability distributions, they can be combined with extreme value 
distributions to estimate the expected number of years in which flooding exceeds a given elevation, inte- 
grated over a given interval of time. Note that this is different from the expected number of flood events in 
a single year; the question here is not, "what is the probability of a flood of at least height X, given the pro- 
jected sea-level change in 2050?" but, "in how many years between 2000 and 2050 do we expect floods 
of at least height X, given the projected pathway of sea level change?" Table 4 shows the expected num- 
ber of years under each RCP with current "1 -in-1 0 year" (1 0% probability per year) and "1 -in-1 00 year" (1 % 
probability per year) flood events for a selection of sites over 2001 -2030, 2001 -2050, and 2001 -2100. 
Figure 7 shows the expected fraction of years with at least one event at the New York City, Key West, Cux- 
haven, and Kushimoto tide gauges for a range of heights and the same periods of time under RCP 8.5; 
additional tide gauges and RCPs are shown in Supporting Information Figure S10. 

At seven of the nine sites considered (New York, Sewell's Point, Key West, Galveston, San Francisco, 
Kushimoto, and Valparaiso, though not Cuxhaven or Stockholm), the expected number of years with 
current 1 -in-10 year flood events, integrated over the 21st century, is under all RCPs at least five times 
larger than the 1 0 that would be predicted without sea-level rise. At the same seven sites, the expected 
number of years in the 21st century with current 1 -in-1 00 year flood events is at least four times higher 
under RCP 2.6 and at least 8 times higher under RCP 8.5 than the 1 that would be expected without 
sea-level rise. 

The increase in expected flood events is influenced both by the magnitude of projected LSL rise and 
by the range of past flood events. The latter is reflected in the difference between the 1 -in-1 0 year 
and 1 -in-1 00 year flood elevations, which will be larger at tide gauges that have experienced more 
extreme flood events. New York City and Cuxhaven are projected to experience fairly high sea-level rise 
(likely 0.7- 1.3 m and 0.6-1 .1 m by 2100 under RCP 8.5, respectively) but have also historically experi- 
enced large flood events, with the 1 -in-1 00 year flood level about 70 cm higher than the 1 -in-1 0 year 
flood level. Under RCP 8.5, these two sites respectively expect nine and four 1 -in-1 00 year floods over 
the 21st century — the same as would be expected for 1 -in-1 1 year and 1 -in-25 year events without 
sea-level rise. 

Stockholm has experienced fairly few large flood events, with the 1 -in-1 00 year flood level only about 20 
cm higher than the 1 -in-10 year flood level, but also has a low projected sea-level rise (likely -0.2 to 0.5 m). 
As a consequence, it also expects nine 1 -in-1 00 year floods over the 21 st century under RCP 8.5. Key West, 
by contrast, has a projected sea-level rise similar to Cuxhaven but has not experienced as many large flood 
events. The 1 -in-1 00 year flood level there is only about 23 cm higher than the 1 -in-1 0 year flood level. 
Accordingly, it is expected to experience 48 years over the 21 st century with a 1 -in-1 00 year flood event, 
about the same as would be expected for a 1 -in-2 year event without sea-level rise. 

The most extreme case among the nine sites considered is Kushimoto, which both has a large projected 
sea-level rise (likely 0.8-1 .3 m by 2100) and has experienced few large flood events, with the 1 -in-1 00 year 
flood level just 10 cm higher than the 1 -in-10 year flood level. Over the course of the 21st century, under 
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Table 4. Expected Number of Years with Flood Events of a Given Height Under Different RCPs a 


1 -in-10 Year Events 


1 -in- 100 Year Events 

Height (m) 

No SLR RCP2.6 RCP4.5 RCP8.5 

Height (in) 

No SLR RCP2.6 RCP4.5 RCP 8.5 


2001-2030 


New York 

1.11 

3.0 

4.8 

4.7 

4.9 

1.80 

0.3 

0.4 

0.4 

0.4 

Sewell's Point 

1.1? 

3.0 

6.0 

5.9 

6.0 

1.66 

0.3 

0.6 

0.6 

0.6 

Key West 

0.43 

3.0 

11.7 

11.5 

11.7 

0.66 

0.3 

0.7 

0.7 

0.7 

Galveston 

0.99 

3.0 

4.9 

4.8 

4.9 

1.89 

0.3 

0.4 

0.4 

0.4 

San Francisco 

0.67 

3.0 

5.8 

5.7 

5.7 

0.88 

0.3 

0.6 

0.6 

0.6 

Cuxhaven 

4.14 

3.0 

3.5 

3.5 

3.5 

4.85 

0.3 

0.4 

0.4 

0.4 

Stockholm 

0.81 

3.0 

2.7 

2.5 

2.6 

1.03 

0.3 

0.3 

0.3 

0.3 

Kushimoto 

1.24 

3.0 

10.7 

10.5 

10.7 

1.34 

0.3 

3.0 

2.9 

3.1 

Valparaiso 

1.17 

3.0 

8.3 

8.2 

8.4 

1.24 

0.3 

1.8 

1.8 

1.9 


2001-2050 


New York 

111 

5.0 

11.8 

11.8 

12.7 

1.80 

0.5 

0.9 

0.9 

0.9 

Sewell's Point 

1.12 

5.0 

15.5 

15.6 

16.1 

1.66 

0.5 

1.5 

1.5 

1.5 

Key West 

0.43 

5.0 

31.1 

30.7 

31.0 

0.66 

0.5 

3.4 

3.6 

4.3 

Galveston 

0.99 

5.0 

14.4 

14.2 

14.9 

1.89 

0.5 

0.8 

0.8 

0.8 

San Francisco 

0.67 

5.0 

17.5 

17.4 

17.8 

0.88 

0.5 

2.1 

2.2 

2.3 

Cuxhaven 

4.14 

5.0 

6.9 

6.9 

7.1 

4.85 

0.5 

0.8 

0.8 

0.8 

Stockholm 

0.81 

5.0 

5.6 

4.5 

5.6 

1.03 

0.5 

0.8 

0.5 

0.8 

Kushimoto 

1.24 

5.0 

79.2 

?9.0 

29.4 

1.34 

0.5 

14.8 

14.8 

15.8 

Valparaiso 

1.17 

5.0 

22.1 

22.0 

23.3 

1.24 

0.5 

7.6 

7.9 

8.9 


2001-2100 


New York 

1.11 

10 

50 

53 

56 

1.80 

1 

4 

6 

9 

Sewell's Point 

1.12 

10 

61 

62 

64 

1 .66 

1 

11 

14 

19 

Key West 

0.43 

10 

81 

81 

81 

0.66 

1 

40 

43 

48 

Galveston 

0.99 

10 

62 

62 

64 

1.89 

1 

4 

5 

8 

San Francisco 

0.67 

10 

65 

66 

67 

0.88 

1 

26 

31 

36 

Cuxhaven 

4.14 

10 

21 

22 

27 

4.85 

1 

3 

3 

4 

Stockholm 

0.81 

10 

15 

13 

23 

1.03 

1 

4 

2 

9 

Kushimoto 

1.24 

10 

79 

79 

79 

1.34 

1 

63 

64 

65 

Valparaiso 

1.17 

10 

68 

69 

72 

1.24 

1 

45 

48 

54 


a Heights for U.S. sites are with respect to the local mean higher high water datum for the 1983-2001 epoch. 
Heights for non-U. S. sites are with respect to the local mean sea level datum for the 1 983-2001 epoch. 


all RCPs, Kushimoto is expected to experience more than 60 years with flooding exceeding the current 
1 -in-1 00 year flood level. 

Sea-level rise allowances [Hunter et a!., 201 3] quantify the amount by which a structure needs to be raised 
so that its current flood probability remains unchanged. For example, the U.S. National Flood Insurance 
Program's Special Flood Hazard Areas are defined as areas with a 1% per year flood probability [National 
Flood Insurance Program, 201 3]. A corresponding sea-level rise allowance would indicate the height above 
the current 1 -in-1 00 year flood zone that would maintain an average 1 % per year flood probability over 
the period of interest. Note that, because the allowance is with respect to flood risk integrated over time, 
its magnitude is less than that of the sea-level change expected by the end of the period of interest. At 
New York City, a project with a 2001 -2030 lifetime, such as a house with a 30-year mortgage, would need 
to be elevated by 1 7 cm above the no-sea-level-rise 1 -in-1 00 year flood zone to maintain a 1 % per year 
flood probability. An infrastructure project with a 2001 -2050 lifetime would need to be raised 26 cm, 
while a project with a 2001 -2100 life time would need to be elevated by 52-69 cm, depending on the 
emissions trajectory (Figure 7). 
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New York Key West 




(c) Cuxhaven (d) Kushimoto 




Figure 7. Expected fraction of years with flooding at tide gauges in excess of a given height under stationary sea level (black) and 
RCP 8.5 over 2001-2030 (blue), 2050 (green) and 2100 (red). Grey vertical lines indicate the current 1 -in-10 and 1 -in-1 00 year flood 
levels. Heights are relative to mean higher high water for U.S. sites and mean sea level for non-U.S. sites. 


6. Cautions 

In addition to highlighting the sensitivities and research needs noted in section 4, we raise several cau- 
tions in interpreting our projections. 

First, in the near-term, internal variability in sea-level rise [e.g., Bromirski etoi, 201 1] makes estimation 
of precise timing of LSL change difficult. Most sites experience interannual variability with a 2 n range 
of about 4- 10 cm [Hay et at., 2013; Kopp eta!., 2010]. At the illustrative sites we consider, the difference 
between the 1 7th and 83rd percentile projections exceeds the decimeter level between 2030 and 2050. 
Until this threshold is reached, year-to-year variability will be comparable to the uncertainty in projections. 

Second, as previously noted, historically estimated background rates of local, non-climatic processes may 
not continue unchanged. For example, while we project 72 ±5 cm of 21st century sea-level rise due to 
non-climatic factors at Grand Isle, Louisiana, changes in fluid withdrawal could reduce the projection 
[ Blum and Roberts, 2012]. 

Third, our background rate estimates are the result of an algorithm applied to a global database of 
tide-gauge data, with different sites having been subjected to different degrees of quality control. Some 
tide-gauge sites may have experienced datum shifts or other local sources of errors not identified by the 
analysis. We recommend that users of projections for practical applications in specific regions scrutinize 
local tide-gauge records for such effects. 

Fourth, our flood probability estimates should be viewed indicatively. They are based on hourly 
tide-gauge records that may be of insufficient length to capture accurately the statistics of rare flood 
events. They do not account for projected changes in tropical or extratropical cyclone climatology, such 
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as the expectation that Category 4 and 5 hurricanes may become more frequent in the North Atlantic 
[e.g., Bender et at., 2010] and perhaps globally [Emanuel, 201 3]. They are developed for specific tide-gauge 
locations where flood risk is likely indicative of, but not identical to, risk for the wider vicinity, due to 
variation in local topography and hydrodynamics. Nonetheless, they do highlight the inadequacy of 
flood risk assessments based on historic flood probabilities for guiding long-term decisions in the face of 
ongoing sea-level rise. 

7. Conclusions 

Assessments of climate change risk, whether in the context of evaluation of economic costs or the plan- 
ning of resilient coastal communities and ecosystem reserves, require projections of sea-level changes 
that characterize not just likely sea-level changes but also tail risk. Moreover, these projections must esti- 
mate sea-level change at specific locations, not just at the global mean. They must also cover a range of 
timescales relevant for planning purposes, from the 30-year time scale of a typical U.S. mortgage, to the 
> 50 year lifetime of long-lived infrastructure projects, to the > 1 century lifetime of the development 
effects of infrastructure investments. In this article, we synthesize several lines of information, including 
model projections, formal expert elicitation, and expert assessment as embodied in the Intergovern- 
mental Panel on Climate Change's Fifth Assessment Report, to generate projections that fulfill all three 
desiderata. 

Under RCP 8.5, we project a very likely mean global sea-level rise of 0.5 -1.2 m by 2100 and 1. 0-3.7 m 
by 2200, which under the strong emissions mitigation of RCP 2.6 is lowered to 0.3-0.8 m by 2100 
and 0.3-2.4m by 2200. Local sea-level rise projections differ from the global mean due to differing 
background rates of non-climatic sea-level change, spatially variable responses to different land ice 
reservoirs due to static-equilibrium effects, and spatially variable ocean steric and dynamic changes. 
Static- equilibrium effects lead to a tendency for greater-than-global sea-level rise in the central and 
western Pacific Ocean. Mid-latitude and high-latitude sites in North America and Europe are generally less 
exposed to climatically driven sea-level change, with the exception of northeastern North America, which 
has potential for a high oceanographic sea-level contribution. At most sites, by the end of the century, 
uncertainty is due primarily to uncertainty in AIS mass loss, though oceanographic uncertainty is also a 
major term at sites where oceanographic processes may make a significant contribution to sea-level rise. 

Probabilistic projections of future local sea-level rise pathways can be combined with statistical or hydro- 
dynamic flood projections to estimate flood probabilities that more accurately assess the risks relevant to 
structures and populations. Projected sea-level rise can dramatically change estimated risks; at the Bat- 
tery in New York City, for example, we project over the 21st century an expected nine years with "1 -in-1 00 
year"' flood events under RCP 8.5 and four under RCP 2.6. Such projections, especially if improved or aug- 
mented by more detailed storm and flood models that include factors such as changes in tropical and 
extratropical cyclone climatology and by hydrodynamic models of overland flooding, can guide insurance, 
land use planning, and other forms of coastal climate change risk management. 
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Erratum 

In the originally published version of this article, the file names in the supporting information readme 
did not match the file names of the supporting information files. The readme content has since been cor- 
rected and this version may be considered the authoritative version of record. 
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